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In heavy ion collisions resonances can create strong non- 
Gaussian effects in the 2-pion correlation data. Hence, the 
commonly used Gaussian fit parameters do not fully charac- 
terize these correlators. We suggest a different set of HBT 
parameters which does not presuppose a particular shape of 
\^ie correlator and allows to extract additional (non- Gaussian) 
C&iformation. Within a simple model for an expanding source 
C&lcluding resonance decays it is shown that this additional in- 
' tdrmation provides a clean distinction between scenarios with 
+^d without transverse fiow. 
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The aim of Hanbury-Brown/Twiss (HBT) interfero- 
r^etric analyses of heavy ion collisions is to extract from 
IShe measured 2-particle correlators C(q, K) of identi- 
£^1 particles as much information as possible about the 
3a3atio-temporal distribution S{x, K) of the particle emit- 
^«ng sources in the collision region. It is based on the 
r*plation gl 
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far, the interplay between the experimentally measur- 
rable momentum correlator C(q, K) and the theoretical 
^ ^ncept of space-time emission function S{x, K) was in- 
vestigated mainly in the context of Gaussian approxima- 
yi(ons for both. Experimental data are commonly fit to a 
. .Gaussian ansatz iSl 
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i,j = o,s,l. (2b) 



Here, the relative momentum components are defined 
parallel to the beam (I = longitudinal), parallel to the 
transverse component of K (o = out), and in the re- 
maining third direction (s — side). Ros = Rsi = for 
azimuthally symmetric collisions. 

The fit parameters Rfj(K) can be calculated from a 
Gaussian ansatz for the emission function S{x, K) in 
terms of space-time variances 



S{x,K) ~ N{K)S{x,K)e'^ 
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In the resulting expression for the correlator, 

C(q, K) = 1 + exp {~q^q''(x^i,)) , (4) 

only three of the four relative momentum components 
are independent, the fourth being fixed via the on-shell 
constraint = q • (K/i^") = q • /3. In the Gaussian 
framework, the consequences of the on-shell constraint 
are well studied. Especially, the HBT radius parameters 
read 



4.(K) = ((i,-/3.i)(i,-/3,f)) 



(5) 



This is the starting point of most spatio-temporal inter- 
pretations of correlation data. Note however that (^ 
presupposes a Gaussian shape of C(q, K) since (^) de- 
termines the curvature components of the correlator at 
q = §: 
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and these coincide with the experimentally determined 
half widths of C(q, K) only for Gaussian correlators. 
Non-Gaussian characteristics of the correlator are diffi- 
cult to quantify and control in this Gaussian framework. 

Here, we suggest to characterize the correlator with 
a different set of a few HBT parameters, so-called q- 
moments, whose extraction does not presuppose a par- 
ticular shape of C(q, K). 

The HBT parameters used in the Gaussian framework 
are based on space-time variances (x^x^), i.e. on ex- 
pectation values {f(x)) of the emission function in co- 
ordinate space, cf. (pq). In contrast, we suggest here to 
base a quantitative analysis of C (q, K) on q- variances 
which are expectation values ((.g(q))) of the true correla- 
tor C(q, K) — 1 in relative momentum space. 

For a Gaussian correlator, the HBT parameters can be 
determined either by fitting to the Gaussian (j^), or by 
computing the integrals 



Jd^qq^qj [C(q, K) - 1] 
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For a non-Gaussian correlator, we can define the HBT 
radius parameters and the intercept A in terms of these 
"(7- variances" . 

The deviations of the correlator from a Gaussian shape 
( "non-Gaussicities" ) are then quantified by higher order 
q-moments. Due to the q— > — q symmetry of C(q, K) the 
first non-Gaussian contribution shows up in the fourth or- 
der q- moments ("kurtosis"). All higher order g-moments 
can be calculated as derivatives of the generating func- 
tion Z(y,K), 



Z(y,K)= J Sqe^^y [C(q, K) - 1] 



lnZ(y,K) 



(9) 
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^From this generating function, the correlator can be re- 
constructed completely. The series of rt-th g-varianccs 
( p^ is merely a convenient way to characterize its shape 
starting with its "Gaussian" widths {{qiqj)) and going for 
increasing n step by step to finer structures. 

To calculate Z(y,K) directly from a given emission 
function it is convenient to use the normalized "relative 
distance distribution" 



p{u; K) = / d'^X s{X + f , K) s{X -^,K), 



d^u p{u; K) = 1 . 
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written in terms of the normalized emission function 
s(x,K) = S{x,K)l J d'^x S{x, K). p is real and even in 
u. Then 

Z(y,K)= Jd^qe'^i y J d'^ue"^ ''piu;K). (12) 

In the Cartesian parametrization, this expression simpli- 
fies to a one-dimensional integral: 



(13) 



Z(y,K)= J dtp{ys,yo + P±t,yi+l3it,t;K). 



For Gaussian emission functions Eqs. (0), (|l^ and (|lJ 
reproduce the relations (||). 

The method of g-variances generally requires to in- 
vert the matrix (|^) and to discuss a 4-dimensional tensor 
of fourth order moments. Such a full analysis will be 
published elsewhere 0. Here, we restrict ourselves to 
a unidirectional analysis of the correlations C((7i,K) = 
C{qi,qj^i=0,K) along the three Cartesian axes. For 
the relative momentum component qi , the corresponding 
HBT radius parameter and intercept are then defined via 
the relations (we use the same notation as for the three- 
dimensional g-moments) 



i?f (K) = 
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Jdq,[C{q,,K)-l] 
A,(K) = (i?,(K)/V^) / dq, [C{q,, K) - 1] . (14c) 



Even Gaussian correlators have non-vanishing higher q- 
moments: 
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Being expressible in terms of the second moments they 
do not contain new information. The non-trivial higher 
order information is contained in the (normalized) q- 
cumulants, in which these trivial contributions are sub- 
tracted: 
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The normalization removes the dependence on the Gaus- 
sian widths of C(q, K) and turns them into dimensionless 
measures for deviations from a Gaussian shape. 

To extract the moments ((gf)) from data one replaces 
(14b) by a ratio of sums over bins in the gj-direction. The 
higher the order n of the g- moment, the more sensitive 
are the extracted values to statistical and systematic un- 
certainties in the region of large qi. First investigations 
with event samples generated by the VENUS event gen- 
erator indicate that the current precision of the data in 
the Pb-beam experiments at the CERN SPS permits to 
determine the second and fourth order g-moments 
Accordingly, we restrict our discussion of non-Gaussian 
features to the "kurtosis" 
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We have computed the unidirectional Gaussian param- 
eters Ri and the kurtosis for a model pion emission 
function including resonance decay channels R: 

S^{x,p) = Sf'{x,p) + J2 Sr^Ax,p) . (18) 



The contributions Sr^^^ are obtained from the direct res- 
onance emission function S'^^'^{X,P) by propagating the 
resonances of widths F, produced at (X^,P^), along a 
classical path x'' = X^^ + according to an exponen- 
tial decay law |H,||: 

Sr^^{x-p)^M ^ dsg{s) J ^SiP-p-ME*) 

xJd^X J drFe-r-5(4)(^_ (x + ^r)) 

xSf{X,P). (19) 

We consider isotropic decays in the resonance rest frame 
with decay phase space g{s), E* being the energy of the 
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observed decay pion in this frame and s the squared in- 
variant mass of the (n — 1) unobserved decay products. 

Our model assumes local thermalization at freeze-out 
and produces hadronic resonances by thermal excitation. 
For particle species i with spin degeneracy 2Ji + 1, the 
emission function reads 01 



SfH^,P) ^^^P-n{^) cxp H{x) , 



H{x) = cxp — 



2i?2 2(Ar;)2 2(Ar)2 



(20) 



The Boltzmann factor exp[— (P-u(a;)— /ii)/r] implements 
both the assumption of thermalization, with temperature 
T and chemical potential fii, and collective expansion 
with hydrodynamic flow 4- velocity M^(a;). Space-time is 
parametrized via longitudinal proper time t = y/f^ — z^, 
space-time rapidity ?/ = ^ In \(t z)/{t ~ z)], transverse 
radius r and azimuthal angle (j). The Gaussian factors in 
H{x) specify the spatial extension of the source as well 
as a Gaussian average around a mean freeze-out proper 
time To with dispersion At. For freeze-out along proper 
time hyperbolas P-n(x) ~ M± cosh(F — ??) @ where Y 
and M± are the rapidity and transverse mass associated 
with P. 
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FIG. 1. Two-pion correlations for transverse flow rjf — 0. 
Curves are calculated for the model emission function 
without resonance decay contributions (dashed lines), includ- 
ing pions from the shortlived resonances p, A, K* , E* and 
from the uj (dotted lines), and including also pions from the 
longlived resonances r], rj' , Kg, E, A (solid lines). 

In the longitudinal direction we assume scaling expan- 
sion, vi = z/t OT r]i — ^ ln[(l -I- vi)/{l — vi)] = rj. For the 
transverse expansion we choose a linear rapidity profile: 



Vt{r) = Vf 



(21) 



All numerical results presented here are obtained for 
the set of source parameters T — 150 MeV, R = 5 fm, 
A77 = 1.2, To = 5 fm/c. At = 1 fm/c and piB = /is = 0. 
We include all pion decay channels of p, A, K*, E*, uj, 



r], r]', K'LYI and A with branching ratios larger than 5 
percent |lO| . 

For this model we have numerically computed the 2- 
pion correlator. Typical results are shown in Fig. |l|. A de- 
tailed model study of how resonance contributions affect 
the shape and pair momentum dependence of the correla- 
tor will be published elsewhere Here, we merely ob- 
serve that due to their exponential decay law, resonance 
decays can contribute exponential tails to the emission 
function Sr^t^ in (|l9|), thereby leading to a generically 
non-Gaussian shape of the correlator. 

To illustrate the use of q-variances we now discuss 
how the space-time characteristics of the emission region, 
specified in our model by the geometric input parame- 
ters and the transverse flow 77/, show up in these new 
2-particle observables. First, we checked that for the 
model (p^pO|) the inverted second g-moments Ri coin- 
cide in all three directions i — 0, s,l with the Gaussian 
fit parameters extracted for a set of n equidistant points 
qp' between and 50 MeV by minimizing 



^(lnC(gp\K)-logA + i?2g^ 



■2^0)' 



(22) 



These "Gaussian" radius parameters Ri and their trans- 
verse momentum dependence is well-studied for model 
emission functions ( po| ) not including resonance decay 
contributions ||^. Especially, it is known that the trans- 
verse homogeneity length Rt of ([20|), i.e. the size of the 
effective emission region in the transverse plane, shrinks 
with transverse flow r/f approximately like |^ 



T 



(23) 



This introduces an ryj-dependent X_L-slope of Rs which 
was interpreted previously as a signature of transverse 
flow. Fig. H shows that if resonance decay contributions 
are added in the calculation, this distinctive behaviour is 
lost: Rs develops a K±-s\ope even for scenarios without 
transverse flow. The reason is that resonances can prop- 
agate outside of the thermally equilibrated region before 
decaying. This effect ("lifetime effect") leads to expo- 
nential tails in Sr^tt and tends to increase the Gaussian 
widths Ri . Since the relative abundance of resonances is 
larger for small K±, this tail is more pronounced in the 
region of small K±, thereby leading to a if^-slope of Rs- 
In contrast, for a finite transverse flow rjj = 0.3 the 
if^-dependence of Rs is due to the flow, and resonance 
decay contributions do not change the slope of Rs in our 
model. This can be traced back to the shrinking trans- 
verse size (2^ of the direct resonance emission function 
S"^"" in (^): due to their larger rest mass, parent reso- 
nances have a smaller effective emission region than ther- 
mal pions in the transverse direction. For the case de- 
picted in Fig. 1^, this effect counterbalances the lifetime 
effect almost exactly. 
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According to (|23|), the absolute size of the side radius 
depends not only on ry/ but also on the input parameter 
R for which no independent measurement exists. As a 
consequence, the Gaussian radius parameter Rs{K±) by 
itself does not allow to distinguish scenarios with and 
without transverse flow once resonance decays are in- 
cluded. This conclusion extends to the other Gaussian 
widths, Ro and i?/, as well. 
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cial additional information on the correlator which can 
help to distinguish physical scenarios which are difficult 
to disentangle on the level of "Gaussian" HBT radius pa- 
rameters. On the other hand, they can still be directly 
related to specific features of the emission function and 
thus should aid the spatio-temporal interpretation of cor- 
relation measurements. This is less straightforward if one 
only knows the overall shape of the correlator, without 
decomposing it into its moments. 

In summary, we expect these new HBT parameters, 
both the one- and multi-dimensional g-moments of the 
correlator, to become valuable additional observables in 
future comparisons between models and experiment. 

This work was supported by BMBF, DFG and GSI. 
Discussions with P. Foka, H. Kalechofsky, B. Lasiuk, 
M. Martin, H.-P. Naef, L. Rosselet, P. Seyboth and S. 
Voloshin are gratefully acknowledged. 



FIG. 2. The inverted q-vanance Rs of (2.8a) and the kur- 
tosis As of (2.11) in the side direction at Y — as functions 
of K±. Left: r]f — (no transverse flow). Right: r]f = 0.3. 
The difference between the dashed and solid curves is entirely 
dominated by lu decays. 

However, the physical origin of the if j_ -slope of Rs 
is different for the two situations shown in Fig. ^. This 
shows up in the kurtosis of the corresponding correlators: 
without transverse ffow, resonance decay contributions 
increase Rs due to the lifetime effect which generates ex- 
ponential tails in Sr^tt- These tails manifest themselves 
in the correlator through a significant non-Gaussian com- 
ponent. For non-zero transverse flow, on the other hand, 
the i^x-slope arises from the ifj_-dependent shrinking 
( p3| ) of the effective transverse emission region. This ef- 
fect is more prominent for resonances than for thermal 
pions, i.e. 5'*'' is spatially more extended in the trans- 
verse plane than Sfl^ , and hence it "covers" a substantial 
part of the exponential tails of S r^tt . As a consequence, 
the total emission function ( p^ ) can be expected to show 
much smaller deviations from a Gaussian shape for the 
scenario with transverse flow, and this again leads to a 
more Gaussian correlator. 

This explains why the kurtosis As{K±) plotted in 
Fig. H provides a clearcut distinction between the two 
scenarios. Aj simply reflects the importance of resonance 
decays. For vanishing transverse flow. As is significant 
and decreases with increasing K± since the relative abun- 
dance of resonance decay pions dies out; for finite trans- 
verse fiow 77 J = 0.3, it is an order of magnitude smaller. 
A iiT^-dependence of Rs without measurable kurtosis is 
thus a clear sign of transverse fiow. 

Here, we do not discuss to what extent this particular 
7;/-dependence of the kurtosis As is generic for realistic 
models of heavy ion collisions including resonance decays. 
This will require a more systematic study p^ ] . We simply 
conclude that higher order g-variances can provide cru- 
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